Localization of spinons in random Majumdar- Ghosh chains 
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We study the effect of disorder on frustrated dimerized spin-1/2 chains at the Majumdar-Ghosh point. The 
Majumdar-Ghosh condition is first generalized in the presence of random bonds, defining a random Majumdar- 
Ghosh (RMG) point with locally correlated disorder. Using variational methods and density-matrix renormal- 
ization group approaches, we show that spinons, the deconfined fractional elementary excitations of the chains, 
remain gapped at the RMG point but get localized in Lifshitz states whose localization length is analytically 
obtained. Increasing the disorder strength at the RMG point induces a transition to a paramagnetic phase of lo- 
calized spinons. Moving away from the RMG point by uncorrelating the random bonds, an effective interaction 
between spinons brings the chain into a gapless and partially polarized phase for arbitrarily small disorder. 



PACS numbers: 75.10.Kt, 75.40.Mg, 75.10.Jm, 75.10.Pq 
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Spinons are fractional excitations corresponding to half of 
a spin excitation in quantum magnets. They typically appear 
in understanding the excitation spectrum of one-dimensional 
systems such as the frustrated J\ — J2 Heisenberg chain. 
This model possesses an exact ground-state at the Majumdar- 
Ghosh (MG) point Ji = 2 J2 [1] which is the prototype of 
a valence bond solid (VBS) state and for which a variational 
approach describes well elementary excitations J2]. Further, 
spinons play a crucial role in unconventional two-dimensional 
phase transitions in which they could be deconfined fH- In- 
vestigating the effect of disorder on their dynamics is essen- 
tial, all the more since randomness is inherent to experimental 
samples. Possible strategies to study random quantum mag- 
nets are bosonization [4], provided the clean system is gapless 
and the disorder is small, or real-space renormalization group 
(RSRG) [5], suited for the strong disorder regime. The lat- 
ter is asymptotically exact in the case of an infinite-disorder 
fixed point [6], but its outcome can be questioned at small dis- 
order when it converges to a finite-disorder fixed point. Nu- 
merical approaches are challenging due to strong finite-size 
effects from rare events |01 ■ Unfortunately, the interplay be- 
tween frustration and disorder cannot be addressed using the 
powerful quantum Monte-Carlo method because of the sign 
problem. Lastly, most studies on random magnets focus on the 
ground-state while little is known about the fate of elementary 
excitations. So far, the ground-state physics of random frus- 
trated dimerized chains has only been addressed within the 
RSRG picture |8[] and has been found to belong to the class of 
the large-spin phase J9|]. 

In this Letter, we unveil some remarkable effects of ran- 
domness on frustrated dimerized chains using a variational 
approach supported by density-matrix renormalization group 
(DMRG) calculations lioll . They provide both quantitative 
predictions and an intuitive picture of the physics. The anal- 
ysis follows two steps. First, the Majumdar-Ghosh condition 
is generalized to random couplings leading to a locally cor- 
related disorder condition dubbed as the random Majumdar- 
Ghosh (RMG) point. There, an Anderson model governs the 
dynamics of a spinon which leads to its localization. Increas- 
ing disorder, a transition to a paramagnetic phase of localized 



spinons is found. Moving away from the RMG point by un- 
correlating the couplings generates an effective interaction be- 
tween spinons which stabilizes the formation of domains and 
leads to the destruction of the spin gap. 

The RMG point - We consider a frustrated dimerized spin- 
1/2 chain with random nearest-neighbor couplings ai and 
next-nearest neighbor couplings 



(1) 



where S$ are spin-1/2 operators. In the following, the average 
couplings are written al — a and = (3 with a = 2(3 to start 
from the MG point. Applying % on the MG state |MG) = 
|-« •-) withdimers \—) = 775 [|t4) - lit)] 

starting on even sites 2j gives 

h\mg) = \J2 (~ 3c ^ + a2 J+ 1 - 0« - hj+x) l MG > 

j 

+ ( a2 ^! - hi ~ hj-i) |2j - 2, 2j + 1) , 

j 

where we write \2j — 2, 2j + 1) = |-« *** ^ •** •-■), 

the state with a dimer on bond (2j — 2, 2j + 1). A similar ex- 
pression is obtained for the other MG state. Thus, MG states 
remain degenerate eigenstates of the Hamiltonian provided 



cti — ^ + f3 i+ i (RMG condition) , 



(2) 



which is a generalization of the MG point in the random case. 
It imposes a local correlation between the random couplings. 
To understand the nature of the ground-state of Hamiltonian 
©, we first investigate the elementary spinon excitations. 

Localization at the RMG point - An effective model for the 
dynamics of a spinon is obtained assuming (0 by considering 
a chain with an odd number of sites and spinon states \j) = 
|-» •** t •** •■■) with a free spin at site 2j + 1 separating two 
MG domains. Following the variational approach [11] which 
consists of projecting H orthogonally onto the free family of 
non-orthogonal states {\j)}, the spinon Hamiltonian reads 



H 



spinon 



li) 



|j-i> + di> + |j + i) 



(3) 
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FIG. 1. (color online) At the RMG point : (a) Lyapunov exponent 7 vs energy E of the effective model, (b) magnetization profile from DMRG 
in sector S z = 1 /2. Inset: log plot, (c) dimerization profile showing the MG domains, (d) finite size effects on the Lyapunov exponent. 



in which the MG state energy £mg 



J2i Pi has been sub- 



tracted. Consequently, the motion of the spinon, taking place 
either on odd or even site sublattices, obeys a special kind 
of Anderson Hamiltonian expected to induce localization, 
^spinon is non-hermitian because {\j)} is not orthogonal but 
this matrix is similar to an hermitian matrix fiUl . The spinon 
wave-function can be written |V'spinon) = Y^j V'jli)' w i*h co- 
efficients ipj. Following the Dyson-Schmidt method |l2, 131. 



we introduce = p2j+i and the Riccati variables Rj = 



tpj+i/tpj to rewrite Schrodinger's equation 

Rj + 5/2 + l/Rj-! = 2E/f3 j 



(4) 



for a given spinon energy E. The integrated density of states 
N(E) and Lyapunov exponent j(E) of a single spinon ex- 
citation are obtained by extending the energy to the complex 
plane. Assuming that the probability density of the Riccati 
variables converges toward an invariant distribution of mea- 
sure dW(R) as j — > 00, the characteristic function 



Q(z) = J dW{R)\nR 



(5) 



is such that Vt(E + i0+) = j(E) + wr(l - N(E)). These 
quantities can be obtained either numerically or analytically 
from a weak-disorder expansion UM as described below. 
In the non-disordered case, the spinon dispersion relation is 
E{k) = /3(5/4 + cos k) with k the momentum. By introduc- 
ing the variables 

-2ifc%-l-ft(l-%-i) 



Vi 



9j 



(5/4- 



- cos k) (f3j 



Vj- 



l) 



ifij sin k 



(6) 



(7) 



it can be shown that the first term in the expansion in the first 
moment of the g-distribution gives £1 ~ ik — g 2 /2. Special- 
izing to the case of a uniform distribution of width 2A, the 
explicit calculation for energies E<f3/A gives 



j(E) = arcosh 



E 
1 



1 



E 2 



*{E-\p)*-p 



(8) 



The result is compared to numerics on Fig QI a )- In partic- 
ular, we obtain that the spinon localization length £ sp i n on = 
l/'y(E m i n ) in the state with the lowest energy E m [ n = /3 m i n /4 
(where /3 m j n = min/3j), scales as: 



a spinon 



(9) 



Notice that f sp i non cannot be captured by RSRG and is not re- 
lated to the spin correlation length of a MG state. 

Another important outcome of the effective model is that it 
provides hints on the finite-size effects on the spinon energies, 
with consequences on the spin gap and the localization length. 
The lowest energy spinon states correspond to the regime of 
Lifshitz localization, in the tail of the density of states and 
controlled by a rare-events scenario. Adapting Lifshitz argu- 
ment IU3U15I1 . a region of length £ with many f3j close to /3 m ; n 
- i.e. provided \(3j — /3 m j n | < C(Eg — E m \ n ) for all j in the 
region, with C a constant - has its lowest energy of the order 
of Ei ~ E min + /3 m i a Tr 2 /2£ 2 , assuming /? min > and writing 
E m i n — /3min/4. Since the probability of creating such region 
scales as Pi oc [C(Eg— E m i n )/2A] e forauniform distribution, 
the low-energy behavior of the integrated density of states is 



N(E) oc exp ■ 



JM^J 2A/C 



E — E„ 



E — E„ 



(10) 



This behavior is in very good agreement with the numer- 
ics on the effective model 111 ill . Regarding finite-size ef- 
fects on the spinon energy £i in a chain of length L, the 
probability to have the minimum energy must be such that 
P e ~ 1/L. This yields N(E L ) - 1/L in Eq. OB, again 
in good agreement with numerical results ifTlll . Asymptoti- 
cally, we thus expect finite-size corrections of the form El ~ 
Emm + if/? m in(lri(ln L) j In L) 2 with K a constant. 

In order to validate this effective model, we compare it to 
accurate DMRG calculations of the magnetization profile in a 
chain with total spin S z — 1/2. To do so, one relates the local 
magnetization rrij = (Sf) to the ipj using ni2j+\ — | IV'jl 2 
and m,2j — —2\ipj\ 2 , assuming that the ipj vary smoothly 
enough. As judged by the results of Fig. |TJb), the effective 
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FIG. 2. (Color online) At the RMG point - typical magnetization 
(a) and dimerization (b) profiles in the lowest triplet excited state for 
A < p. (c) evolution of the singlet and triplet spin gaps, as well as 
the spinon density iV(0), v.s. the disorder strength A. 



model provides quantitative predictions of the magnetization 
profile. From the local dimerization pattern <i, = (S^ • Si+i) 
of Fig. QI C X the localized spinon clearly separates two dif- 
ferent MG domains. One can extract the actual £ S pj non from 
DMRG profiles. On Fig. Hid), strong deviations between 
DMRG calculations and the infinite size result (0 are ob- 
served. This difference actually originates from finite-size ef- 
fects : one has to take £ sp i non = \/^{El) on a finite system, 
which yields strong deviations, even for large sizes. 

Spin gap and transition to a paramagnetic state - From 
these results on single spinon excitations, we can infer the 
behavior of the spin gap at the RMG point with increasing 
disorder. The MG states remain degenerate ground-states at 
small enough disorder and the lowest triplet excitation above 
them is to create two localized spinons far away from each 
other with both the minimal energy E m i n - They naturally be- 
long to two different sublattices and are expected not to in- 
teract on sufficiently large chains. Consequently, the triplet 
gap, within the effective model approach, is expected to be 
A 5 = /3min/2 in the thermodynamical limit and degenerate 
with the singlet gap. We plot in Fig. [2] the evolution of the 
DMRG triplet and singlet gaps at the RMG point for a box 
distribution, with increasing A. In the clean case, the vari- 
ational approach slightly overestimates the gap, in particular 
due to the core interaction between deconfined spinon, and 
DMRG finds a gap smaller than (3/2. The effective model 
prediction reads As = (f3 — A)/2 and vanishes for (3 = A, in 
relatively good agreement with the data. We also give A$(L) 
obtained from the effective model to show that DMRG data 
are essentially controlled by the strong finite-size effects dis- 
cussed above. Actually, the fact that spinons are localized 
and thus hardly interact improves the validity of the effective 
model at intermediate disorder. 

Increasing further the disorder strength, E m [ n becomes neg- 
ative (Emm = 9/3 m i n /4 for /3 m ; n < 0) so that states with two 
or more spinons get energetically favored. The MG states are 
still eigenstates but no longer the ground-states. The resulting 
picture is a paramagnetic phase of localized spinons, neglect- 
ing tiny residual magnetic couplings between spins degree of 



freedom, provided the spinon density is small enough. The 
density of spinons in each sublattice is the density of nega- 
tive energy spinon states A^(0). So is the total spinon den- 
sity. Therefore, the magnetic susceptibility must change from 
zero to x ~ -W(O) across the transition. DMRG calcula- 
tions actually confirm this picture, showing that the ground 
and first excited states are nearly degenerate states with lo- 
calized spinons. The order of the quantum phase transition 
from the gapped to the paramagnetic phase naturally depends 
on the disorder distribution. The location of the transition, 
within the effective model picture, is given by Emm = and 
one expects a continuous transition for a continuous disorder 
distribution. Binary disorder, on the contrary, would yield a 
first order transition. 

In order to determine iV(0) for the effective model, one 
must realize that the Lifshitz argument cannot be used for 
spinon energies close to zero. Indeed, when / 3 m m — 0, the 
energy E# of a state in a cluster of size £ no longer depends on 
I, Yet, one remarks that the number of negative energy states 
is actually given by the number of negative f3j. In the case of 
a box distribution, one obtains that for A > f3: 

N(0) = (A - /3)/(2A) , (11) 

and the associated critical exponent of the susceptibility is 
simply one. Eq. (fTTT i is checked numerically on Fig. |2jc) 
within the effective model calculations. Checking this law us- 
ing DMRG is particularly difficult as the spinon density gets 
very small close to the critical point. Of course, the quantita- 
tive predictions from the effective model are to be taken with 
care because assuming independent spinons requires a low 
enough spinon density and that neglecting states with non- 
local dimer is questionable in this regime. In particular, effec- 
tive magnetic couplings between spinons are non-negligible 
at high spinon density and the phase is likely to be partially 
polarized and not strictly speaking paramagnetic. 

Away from the RMG point - It is now instructive to turn 
to the situation where the RMG condition (O is progressively 
weaken by uncorrelating the a,; and fa . This can be done using 
random variables Si, uncorrelated to the fa, and such that 

ai = (1 - X)(fa + fa+i) + X8i , (12) 

where Si = 2/3 and A e [0, 1] is a tuning parameter. The 
correlations ctifa — a-7/3; = (1 — A)er| (er 2 denoting a vari- 
ance) show that the as and /3s get uncorrelated for A = 1. As 
cr 2 = 2(1 — A) 2 cr 2 + A 2 er 2 , we further impose c 2 = 2cr| to 
study the effect of the correlations only by keeping the disor- 
der strength as for the RMG point. Lastly, using this decou- 
pling, the Hamiltonian (Q]) nicely splits into two parts with a 
more transparent physical input : 

H = Krmg + A 7?jSj ■ S l+1 , (13) 

i 

with rji = Si — fa — fa+i for which we have rjl = 0. Wrmg 
is the Hamiltonian at the RMG studied above. The second 
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FIG. 3. (Color online) (a) Magnetization profile of a MG chain with 
uncorrelated couplings displaying localized spinons (b) Dimeriza- 
tion pattern and corresponding average of the random 77-dimerization 
term of Eq. dl4t over MG domains (orange line). 



term corresponds to a random dimerization term with pos- 
itive and negative couplings. At small A, the second term 
acts as a perturbation which locally lifts the degeneracy be- 
tween the two MG states of Wrmg- Clearly, MG states are 
no longer eigenstates since high order terms in perturbation 
theory put weights on states = |-» trS^i t - •-.) 
with a dimer on bond More than dressing the MG 

state, the random dimerization actually destroys the spin gap 
as soon as A 7^ and erg 7^ 0. Indeed, the energy difference 
SEij = (i } j\H\i,j) — E M c between a long dimer state and 
the corresponding MG state reads |ll ] 



SE^ = -(/3 2i -i+l32j)+\(r)2j-i + 



3-1 

J2 ^ 2 »- 



i-mn . (14) 



The first part is reminiscent of the spinon localization mecha- 
nism induced by Hrmg while the second term can be under- 
stood as an effective long range interaction between spinons 
at i and j stemming from the random dimerization. Then, a 
rare-event scenario supports the fact that for a infinitely large 
system, it is always possible to find a region with 77s that do 
not compensate and such that the random interaction scales 
as the region size and is negative. Then, since creating two 
spinons will not cost more than an energy proportional to j3, 
one can make SE^ < and the spinons localize at the edges 
of the region (not necessarily in Lifshitz states) with a typical 
length scale £ S pinon - the region size can be made much larger 
than Spinon- As £ sp inon is found to be little affected by A, we 
expect that the localization phenomenon remains mostly con- 
trolled by "Hrmg- Physically, the random dimerization locks 
MG domains to minimize their energy and spinons are local- 
ized at their interface. Lastly, as the region can be as large 
as possible, the effective magnetic coupling between the two 
spinons is arbitrarily small so that the spin gap vanishes. 

This argument is not a mathematical proof of the nature 
of the small disordered phase. Manipulating superpositions 
of states with many spinons or long dimers to have quan- 
titative results in this disordered system, as we obtained at 
the RMG, becomes difficult. In order to support this picture, 



we performed DMRG calculation for non-zero A. We do ob- 
serve that, on a finite-system, the spin gap is strongly reduced 
by increasing A, or by increasing A at fixed A. We exhibit 
in Fig. [3ja) a sample with A = 1 where four spinons are 
present in the ground-state and for which the triplet gap is zero 
within numerical precision. In a infinite chain, this pattern and 
similar ones will necessarily appear, making the ground-state 
gapless. Ferromagnetic couplings between two neighbouring 
spinons cannot be excluded so that the phase should actually 
be partially magnetized. The averaged 77-dimerization field on 
Fig-Sb) supports the locking mechanism of the MG domains. 
Again, strong finite-size effects prevent us from studying nu- 
merically the A->0 and A — > limits without bias. 

Comparison with the RSRG approach - The argument de- 
veloped above, based on localized spinons, is only valid for 
low enough spinon density, ie. at small disorder. It is inter- 
esting to compare these results with the RSRG method best 
suited to the strong disorder limit. From the RSRG equations 
given in[ljj, we notice that the degeneracy of the MG domains 
at the RMG point translates into an instability in the RSRG 
procedure. As soon as A ^ 0, the gap distribution in the 
system converges toward an invariant power-law distribution 
with a non-universal exponent, caracteristic of Griffith phase 
similar to previous studies [8]. Indeed, due to frustration, the 
RSRG equations can generate a few effective ferromagnetic 
couplings, building a large-spin phase |9[]. This is in agree- 
ment with the low-disorder picture and this partially magne- 
tized phase would extend from infinitely weak to strong dis- 
order continuously. A last remark is that the dimer basis used 
above for variational calculations has a deep connection with 
the RSRG picture that tries to catch the most probable dimer 
configuration from the coupling distribution. 

Conclusion - This work provides quantitative results on 
the interplay between frustration and disorder in random MG 
chains. We identify two mechanisms at play : the localiza- 
tion of spinons through a non-usual Anderson mechanism and 
the locking of MG domains by a random dimerization poten- 
tial. The immediate destruction of the spin gap upon putting 
disorder is to be contrasted with its robustness for the explic- 
itly dimerized or spin-1 chains which have a non-degenerate 
ground-state fl6ll . The presence of degenerate MG states 
makes the system very sensitive to disorder. We expect the 
same phenomenogy to play a role in other random VBS. The 
localization of spinons could also be captured by numerical 
methods working in the spinon basis 11711 . Lastly, this model 
gives an example of spinons or "free spins" generated by ran- 
dom couplings which could be experimentally relevant and 
distinct from those generated by vacancies or adatoms. 
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Supplementary material for: Localization of spinons in random Majumdar-Ghosh 

chains 

Non-crossing dimer basis and MG states 

We gather some useful results on the non-crossing dimer basis used for variational calculations. The basis is a free basis of 
non-orthogonal states. In the MG physics, the states with dominant weights are rather simple as they are essentially states with 
nearest-neighbor dimers, with possibly slightly longer dimers locally. 

Applying a nearest neighbour term on a MG state : 

Si ■ Si_)_i|--» •*• •-■) = ^|'"* •** •** •") + 2^'"* *** •■■) (15) 

Applying a next-nearest neighbour term on a MG state : 

Si • S^_4_2 1 -- * •** •** •** •** •■■) = ^|~" •** •** •** •** •■■) — 7j|™« * S™5 "* •-), (16) 
which can be rewritten in the non-crossing dimer basis, using 

|... M <W> M .-) = |- — ~ — ~ .-) + I- ~ ~ .-.). (17) 

Derivation of the single-spinon effective Hamiltonian 

In order to obtain the effective Hamiltonian for the motion of a single spinon, one has to project the Hamiltonian on the 
variational basis { | ■■• 4 •■■)}. Applying Si ■ Sj+i on a spinon at position i : 

Si ■ S l+1 \- ~. t M ...) = M t ~. ...) + M ^ i ...) (18) 

Using the following relation 

|... ~ ^ ...) = |... M t ~ ...) + |... „ _ t ..,}, (19) 
one can deduce the application of S; • S^+2 in the variational basis : 

Si ■ S i+2 |- — * ~- •■■) = -~|- — * ~- •■■) - i|- ~. ~. * .-) (20) 
Applying S;_i • S^+i on a spinon at position i : 

Si-i ■ Si+ii— — * — -) = \\- ~ * — •■■) + \\- ..) (2i) 

One can notice that this state is orthogonal to the variational basis. As a result it simply disappears within the variational 
approach. 

Finally, applying the Hamiltonian on a spinon at position i and projecting on the variational basis : 

■Hj-i,i,i+i|--» •■■) =^("i ^ ft+i + a.i-1 - •*•*•*• •■■) (22) 

+ -fa - ~* ~- i •■■) + - /3i_i)|- 4 — — •■■) (23) 

Some features of the Majumdar-Ghosh state in the presence of disorder 

We discuss some remarkable features of the MG state |MG) in the presence of disorder. We write a and j3 the average values 
of the random couplings and assume that a — 2(3. From ( TToT ) and using 

/..» M M M M •-]-• rf^?*^ »~+ •-) = — 1/2 (24) 
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and that Y^ji' " ' ) " ) f° r large enough system 



L/2 



(MG|-H|MG) = S MG = -7 £ «2j-i -§Lj8 



(25) 



Notice that the above energy is exact even when (0 is not satisfied, but that the simplification from averaging only comes with 
the thermodynamic limit. In particular, it is remarkable that the energy remains independent of the disorder strength A. The fact 
that the MG state is not an eigenstate in general (when (O is not satisfied) can be captured by calculating the energy dispersion 
of the state: 



•XT 

a 2 MG ee (MG|^ 2 |MG) - E&q -)■ — [a* + 2aj + 2c? + 8/3 2 - 4(a/3 + a/3)] 
On can check on this expression that we do have umg = when (fJJ is fulfilled, as expected for an eigenstate. 



(26) 



Energy of a long dimer state 



We denote \i,j) = ■-• 



•3 • •- ), the singlet product state with a long dimer between sites 2i — 1 and 2j (i ^ j). 

This state can be seen as a singlet state between two localized spinons. So, away from MG line (A > 0), it may have a lower 
energy than the MG state. Let us apply the Hamiltonian on this state : 



. / i-l j-2 \ 

H\i,j) =~ I -3 ^ a 2 n-l + ^2{a2n - fan - fan+l) ] \h 
\ n=l n=l ) 



1 2 
"2 / ](&2n - fan - fan+l)\-» •- ■- 

n=l 

1 / j ~ X j_1 \ 

^(-3^a2«+ ^2 (a 2n -i- fon-l- fan)] 

\ n—i n—i+l / 

i J'— 1 

-- \ a 2n-l — fan-1 — fan)\-» ■-• **«^» •- •- 

n—i+1 

-j ( -3 ^ a 2 „-i + ^ (a 2n - /3 2n - A^+i) I 



n=j+l 



n=j+l 



^ ("2n - /?2n — #2n+l)|"« *^ 



n=j+l 

+ ^("2i-2 + Q! 2 i-1 - fai- 2 - fai + a 2 j-l + 0>2j - fai-l - faj + \)\i,j) 



(«2i-2 - Pn-2)\i ~ + ^{a 2 i-i - 02i)\i + 1, j) + (a 2 i-i - &y-i)|i, j - 1) + {a 2j - /3 2 j+i)|i, j + 1) 



+\{fai-i + faj)\i,j) + \ 



2i - 1 



3m- i -• 

Of course, this state is not exactly an eigenstate but using 



(*,i|»±l,j) = (*,j|i,i±l) = -5) 



one can calculate its energy : 



L/2 



a 2 n-i 



(27) 



(28) 
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and compare it with the energy of the MG state : 



- (MG\H\MG) = - ( ^ a 2 „-i - ^ a2n ) = J O^- 1 + + H^V-i + Ys^n-i - 



j'-i 



(29) 



We recall that the 77s variables have a mean value 77 = and variance 

1-A 



0% = 4: 



A 



(30) 



Effective model 



Similarity to a symmetric matrix 

The effective model of Eq. (0) is in a non-hermitian form due to the non-orthogonal nature of the dimer basis. The associated 
matrix simply reads (writing s = 5/2): 



H = ~ 
2 



(sfa fa 

fa sfa fa 

fa sfa fa 



\ 



fa sfa 



\ 



(31) 
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When all the fa are positive numbers, using the similarity transform diagonal matrix D = diag(y / /3i, y/fa, . . .) puts H into the 
following tridiagonal symmetric form 



/ sfa 4fafa 
Vfafa sfa Vfafa 



D^HD 



\ 



Vfafa sfa Jfafa 



Vfafa sfa 



\ 



(32) 



clearly showing that all eigenvalues are real, as for the eigenvectors. 

When some of the (3s are negative, applying the same transform then leads to a complex symmetric matrix, but not to an 
hermitian matrix. A non-diagonal similarity transform is then required. Looking at the 2x2 and 3x3 cases shows that it becomes 
pretty difficult to construct a similarity transform matrix. In addition to showing that the eigenvalues are all real, finding such 
transform would have been very useful to get numerically the spectrum for very large systems. Still, we numerically observe 
and expect that all the eigenvalues and ipj coefficients are real since the effective Hamiltonian stems from an initially hermitian 
Hamiltonian. 



Weak-coupling results 

We gather the weak-coupling results on the Lyapunov exponent of the effective model. The energy is parametrized through 
the variables t or k: 



a , 1 f(j ~ coshi)A 



t „f< E <f a n dE = «| + cos t ): 7= i(<I±^y 

ar-i k 1 /(I + cosht)A 

for£>M md ^ (! + cosll , ): 7 = t __(k__J_ 
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FIG. 4. (Color online) (a) Integrated density of states at low-energies compared with Lifshitz argument od Eq. l llOt . (b) Finite-size effects on 
the spinon energy. 



Numerics on the Lifshitz tail 

We give in Fig. Sfa-b) the comparison between numerical calculations on the effective model and finite-size corrections 
obtained from Lifshitz argument. 



A wrong argument for the susceptibility exponent at the random MG point 

A naive argument for /3 m j n < is the following : the susceptibility \ of the paramagnetic phase should correspond to 
independently filling spinons in single-spinon Lifshitz states up to zero energy. One is tempted to use a Lifshitz formula for 
N(0), which is similar to (ITOb with the changes E# ~ E m i n — (3 m i n ir 2 /2£ 2 and E m i n — 9/3 m i n /4. Using /3 m ; n = A c — A, with 
A c = /3 the critical disorder strength, one gets for the susceptibility exponent = n\/2/3 ~ 1.481 for the uniform distribution. 
Actually, such a prediction is wrong for the reason that the Lifshitz formula does not work close to the E = while it does work 
close to E m in- This behavior is true whatever the smallness of the disorder strength. Instead, we observe that N(E) is linear 
close to E = and that N(0) is rather linear with A — A c . As discussed in the main text, the correct N(0) is obtained by 
coupling the number of negative /3s which leads to 4> = 1 for the uniform distribution. 



On the convergence of DMRG calculations 

DMRG calculations were performed using the finite-size algorithm, targeting one or two states (for instance to determine the 
singlet gap) and keeping typically from 400 to 1000 kept states. As MG are products of dimers, they have a simple matrix- 
product form which makes DMRG pretty efficient. The energies are converged to high precision. In the presence of localized 
spinons and in the S z = sector, the local magnetization should be zero everywhere. Yet, DMRG builds up a variational states 
with non-zero magnetization at the place of localized spinons. Indeed, due to the localization of spinons, triplet and singlet 
states are degenerate within an energy gap that is tiny (we observed gaps below 10 _8 /3) and controlled by the residual magnetic 
couplings between spinons. Thus, the effective couplings between spinons become so tiny that it is extremely hard for DMRG 
to differentiate between the singlet or S z — triplet state and gives a superposition of these states as an output, with a finite 
local magnetization. Still, the spinon localization and magnetization profiles in the S z — 1 are very well converged. 
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RSRG equations for the dimerized chain 

Due to frustration, ferromagnetic couplings can be generated during the RSRG scheme so one has to take into account the 
possibility to generate spins higher than 1/2. The renormalized couplings, which are here written in the general form Jy, depend 
on the spin size Sj. We have the following two equations corresponding to the decimation scheme sketched in Fig. [5] 

Jo 2 



S 



Si 



Jl2 



So 




J03 



FIG. 5. Decimation scheme for the RSRG procedure of Model (TJ. 



if si 7^ S2 or J12 < 0, we take s = si + S2 for Jyi < and s = \s% — sa| for Jyi > and we have 

~ s(s + 1) + si( Sl + 1) - s 2 (s 2 + 1) j , s(s + 1) + s 2 (s 2 + 1) - Si(si + 1) 
J — — — Joi H ;r~ ? — TT\ ^02 



2s(s + l) 



2s(s + l) 



(33) 



if Si = S2 and J12 > 0, we have 



7 t 1 2 / 1 is («^01 - ^02)(^23 - Jl3) 

J03 = J03 + t:Si( s i + 1) ; 



o Jl2 



(34) 



In particular for the first decimations, if a, is the strongest coupling, spins i and i + 1 are decimated and the renormalized 
couplings between remaining spins are : 



■/i-1,1+2 
Ji-2,i+2 
Ji-l,i+3 



(ctj_i - ftXaj+i - _ (A-i + A?7,_i)(/3 i+ 2 + AT? i+ i) 

2ai 2«i 
Pi-i(a i+ i - (ij+i) _ Pi-i{Pi+2 + Ar? i+ i) 

2aj 2a l 
(a»_i - Pi)(3 l+ 2 _ {fi%-\ + Ar?i_i)/3 i+ 2 



i-2,i+3 



2a, 

Pi— lPi+2 

2a, 



2a, 



(35) 
(36) 
(37) 
(38) 



Thus for A = 0, we end up with four degenerated couplings, implicitly reminiscent of the degeneracy of the MG domains at the 
RMG point, which makes the continuation of the RSRG procedure unstable numerically and ill-posed. 

We have studied the behavior of the RSRG equations which lead to a large-spin Griffith phase but the details will be published 
elsewhere. 



2i - L. ' 



